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Abstract. The thermodynamic limit of the internal energy and the entropy of the system of 
quantum interacting particles in random medium is shown to exist under the crucial requirements of 
stability and temperedness of interactions. The energy turns out to be proportional to the number 

■ of particles and/or volume of the system in the thermodynamic limit. The obtained results require 
y—i , very general assumptions on the random one-particle model. The methods are mainly based on 

■ subadditive type inequalities. 

<n: 

^ ■ 1. Introduction 

Since the fundamental work |And58j of P. W. Anderson, the theory of random Schrodinger 
; operators has been an extensively studied field of mathematical physics. The greatest attention has 
r ^ 1 - been paid since now to the one-particle approximation and we do not try to list here even the major 
o f works on this topic. 

There are relatively few papers where finitely many particles are considered. There is a series 
of papers by Michael Aizenman and Simone Warzel that generalize the techniques of fractional 

■ moments method (see, for example, |AW09j ) and another series of articles by Victor Chulaevsky, 
—i. Yuri Suhov and their collaborators that make use of the multiscale analysis [CS09, CBdMSllJ. The 

common point of these works is that they consider the number of particles being fixed and study 
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J> ■ the infinite volume limit for such system. 

The present paper is an attempt to give an insight of what happens if both the number of 
! particles and the volume go to infinity together so that the number of particles per unit volume is 
CN " kept constant. 

The same question has already been addressed by various authors in case of absence of background 
■ potential, i.e., when one-particle propagation is given by pure Laplacian. In this paper, we will 
frequently follow the framework developed by David Ruelle in |Rue99] . though the presence of 
random potential presents certain mathematical difficulty, which we will explain later. We would 
also like to refer to an outstanding article |LL72] of Elliot H. Lieb and J. L. Lebowitz, where 
Coulomb interactions (always in absence of background potential) are treated. 

The idea that the number of particles grows with the volume looks natural in the context of 
condensed matter physics. As a reference real-world example consider a piece of metal or semicon- 
ductor. A bigger piece should contain proportionally more electrons. As macroscopic objects are 
composed of many atoms (Avogadro constant Na ~ 6 x 10 23 mol _1 ), and thus, ions and electrons, 
it turns out that the corresponding mathematical notion is the thermodynamic limit. Its existence 
for thermodynamic quantities, such as internal energy, free energy, calorific capacity, and so on, is 
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the mathematical verification of the fact that these quantities are extensive. The latter is barely 
assumed in physics but actually needs rigorous verification. 

Let us briefly discuss the mathematical objects we study. All the notions will be introduced later 
in full regularity. Let 

H u = -A d + V u (1.1) 

be the random Schrodinger operator that describes a single quantum particle in random environ- 
ment V u . Kinetic part is <i-dimensional Laplacian. One may also consider magnetic Schrodinger 
operator or whatever, provided that a number of basic facts, such as Wegner estimate (see Proposi- 
tion ES]), from the theory of one particle random operators hold true. Actually, the whole ideology 
of this paper is that we take one-particle operators for known and deduce on this base properties 
for multiparticle operators. 

The restriction of to the domain A is denoted by H U (A). For one particle Hamiltonian as in 
( II. ip . we define, with a slight abuse of notation, the n particle operator (restricted in physical space 
to domain A) with pair interactions potential U by 

n 

H U (A, n) = -A nd + V u (j) + U ^ ~ xi )> 
i=i i^j 

where x 1 G A, i — 1, . . . , n, are particles' coordinates. 

Using the notations introduced above, the general question we want to understand is the behavior 
of H u (A,n) in the thermodynamic limit: 

H u {A,n) -)> ?, |A| ->■ oo, n ->> oo, n/\A\ ->■ const. (1.2) 

In this paper, we answer a much more modest question than (II. 2p . Namely, let (B w (A,n) be the 
ground state energy of H u (A,n). In Theorem 13. 5[ we show, in particular, that the ground state 
energy per particle admits the thermodynamic limit: 

<E W (A, n) 

3 lim— - — - — , |A| —7- oo, n — > oo, n/\A\ — > const. (1-3) 

n 

Moreover, the same theorem gives a bit more general result that allows to scale on the eigenenergy 
number in the spectrum. Roughly, the eigenenergy number (counting function) should be of order 
of exponent of the number of particles to ensure the convergence. 

Theorem 13.91 gives the reciprocal result interchanging roles of energy and the counting function 
in the spectrum (the theorem is stated in terms of entropy which is the logarithm of counting 
function) . 

The main tool we use to obtain our results is a modified version of subadditive ergodic theorem 
(see Proposition 14. 6p . For instance, one may show that the ground state energy <E U (A, n) is additive 
with respect to the pair (A, n) up to an error term that can be taken into account. To make use of 
subadditivity we follow the construction of D. Ruelle |Rue99] . Nevertheless, significant modifications 
are made in the proof because of the fact that instead of full translation invariance of free Laplacian, 
we have only the covariance property of the family of random operators. In general, we are only 
able to prove the convergence in L 2 with respect to randomness (see Theorem 13.51 case (nj)). A 
stronger convergence in L 1 and almost surely is established for compactly supported interactions 
(Theorem 13.51 case (Tb]l ). 

In the last part of the present paper, we consider the system of noninteracting fermions in random 
medium. We show that nontrivial effects arise due to Fermi-Dirac statistics even in absence of 
interactions. In particular, we give an exact expression for the limit (II. 3p in terms of the one 
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particle density of states measure (see Theorem 15. 13[) and we find an interesting relation with the 
Fermi energy. 

The rest of the paper is organized as follows. The model of interacting quantum particles in 
random media and the notion of thermodynamic limit are introduced in Section [2 The results 
(mainly on the existence of thermodynamic limit) constitute Section [31 followed by the proofs 
in Section HI In addition, Section H] uncovers some extra properties of the energy density (see 
Subsections 14.31 and 14. 4p . The proofs themselves may be instructive as well. In Section simple 
calculations concerning the thermodynamic limit for vanishing interactions are provided. 

The author is grateful to his thesis advisor, Prof. Frederic Klopp, for proposing the problem, for 
his constant interest and support, as well as many valuable discussions. 

2. Model and Notations 

2.1. Model of Interacting Quantum Particles in Random Media. We consider a system of 
n interacting quantum particles in a random medium. The discrete and continuum cases are treated 
simultaneously and an explicit indication is given if a result is valid only for one setting. In the 
discrete case, the configuration space is given by V = 7L d and for the continuous case by V = M. d . 
In uniform manner, the one-particle Hilbert space is given by 

The n-particle Hilbert space definition depends on the statistics (physical nature of quantum par- 
ticles). The following statistics are considered. 

(a) The Maxwell - Boltzmann statistics. The particles are physically distinguishable and 
no restrictions are imposed on a multiparticle wavef unctions. This model is suitable, in 
particular, for the description of heavy atomic nuclei, i. e., for particles that exhibit classical 
properties. The corresponding Hilbert space is given by 

n 

W = (g)i3 = L 2 (V n ). 

(b) The Bose - Einstein statistics: the particles are bosons. The wavefunction is necessarily 
symmetric with respect to the permutations of coordinates: 

^ = Sym"f, = 4(V"), 

where Sym is the symmetrised tensor product. 

(c) The Fermi - Dirac statistics: they describe fermions. Wavefunctions are restricted to the 
antisymmetric subspace 

n 

ST_= f\Sj = L 2 _(V n ), 
i=i 

where /\ is the external product. 

and S^"i are proper subspaces of fj n . For (j G {0, +, — }, we write P$ to denote the orthogonal 
projector on where (j = stands for the Maxwell - Boltzmann statistics, (j = + for the Bose 
- Einstein statistics and (j = — for the Fermi - Dirac statistics. Obviously, P = P = l^n is the 
trivial projector. 

One particle Hamiltonian is given by 

H U> = H U (1) = -A + V U , 

and acts on Dom(if w ) C .fj, where 
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• A is either discrete or continuous Laplacian, 

• a random potential V w is (at least) Z d -ergodic and satisfies a decorrelation (independence at 
a distance) condition: 

3Rq > 0, such that if dist(A,B) > i? ,then {V^{x)} x& a and {V w {x)} x& b are independent. 

(IAD) 

Remark 2.1. We also take into account the classes of random potentials that have the ergodic group 
reacher than Z d -translations. For instance, everything what follows remains true for the Poisson 
model. 

Notation 2.2. We write Q, P and E for the associated probability space, probability measure and 
expectation respectively. For 7 G Z d , we denote by r 7 the corresponding translations (measure 
preserving transformations) in Q and by T 7 the corresponding unitary transformations (coordinate 
shifts) in L 2 (V). Namely, 

H Tj (u) —T*H U T 7: (2.1) 
where T 7 /(x) = f(x — 7), / £ L 2 (V), x G V. 

By ii/<y we denote a corresponding operator in Sj n that acts only on the i-th particle. More 
precisely, 

H® = Isi ® . . . ® I), ®P W ® 1^ ® . . . ® 1^ . (2.2) 

i — 1 times n — i times 

The n-particle Hamiltonian in random environment V w and with interactions W is given by the 
following self-adjoint operator on $$\ 



H u>i (n) = P„ 



(2.3) 



.1=1 

For each n G N, W n is an interaction potential given by a function of the n particles coordinates 
G V. We refer to the whole collection W = {W n } ne ® as interactions in general. 
Remark also that in this model interactions are deterministic and all particles live in the same 
random background potential V u . 
In (12. 3p . the free part 

n 
i=l 

is called the second quantization of H u in context of the Fock space (see, for example, [BR97] ) . 
Namely, we have to restrict the second quantization of H w to the n-particle subspace of the whole 
Fock space: 

H°^(n) = dT(H^ r 
where dT denotes second quantization procedure. 

Remark 2.3. H®^(n) acts from fjj 1 into itself for any choice of (j, whereas an arbitrary interaction 
potential W n does not necessarily preserve complete (anti) symmetry. That is why the projector 
a-priori acts non trivially in this formula. However, potentials that we consider later are permutation 
symmetric (confer Section [31 property (PI)), so that the projector becomes obsolete in ( 12. 3ft . i.e., 

HUn)=H°Jn) + W n . 
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The Dirichlet and Neumann restrictions of H u ^(n) to a finite box A C V are denoted by H*JA, n), 
where * G {D, N}. if* (A, n) is a self-adjoint operator on ^(A) = Lj(A n ). We omit in 
notations frequently. 

The operator H U $(A, n) has a discrete spectrum. We call the counting function associated to this 
operator 

N Uti (E, A, n) = card{£ fc (A, n, w; ft) < -E}, 
where E^(A, n, oj; (j) are the eigenvalues of H U $(A, n). For the reasons that will become apparent 
later, the entropy is a more convenient quantity: 

S uA {E, A, n) = log jVLj(^, A, n). (2.4) 

Notation 2.4. Sometimes we will drop some (if not all) of the indices and arguments of the counting 
function and the entropy. For example, if we are interested in the dependence on energy, we will 
write just: 

S{E) = S^{E,A,n), N{E) = N U ${E, A, n). 

Remark 2.5. As the counting function M takes its values in NU {0}, the entropy takes its values in 
logNU {-oo}. 

Observation 2.6. For fixed u>, (j, A and n, the entropy E \-¥ S u $(E,A,n) is a non-decreasing 
right- continuous step function. 

The monotonicity of S(E) = S^^E, A, n) allows to define a (quasi-) inverse function E^^A, n, S). 
As jV(-, A, n) is not a local bijection at any point, the inverse function doesn't exist in a canonical 
manner. Our choice of the inverse is the following. 

For S such that e s G N we define 

E Utt (A, n, S) = E cxpS {A, n, u; (j). (2.5) 

The application S (->■ E(S) is a right inverse of the entropy (12 .4p in the following meaning. For 
S G log N one has 

S^(E^(A,n,S),A,n) = S. (2.6) 
Reciprocally, if E ^ Ex{H u ^{A : n)), then 

E Utt (A,n,S Uti (E,A,n)) = E-, (2.7) 

where E~ is the closest from below to E eigenenergy of H^^E, A, n). 
The relations ( 12. 6 p and (12. 7p motivate this choice of an inverse function. 

Definition 2.7. We denote by (£ = € w j(A,n) the ground state energy of the operator H u $(A,n): 

e „ j( A, n)= inf (3d£SML . 

<peDom(H u>i (A,n)) \\(p\\ 2 

Two characterizations of the ground state energy in terms of entropy are given below. 

Proposition 2.8. € is the ground state energy if and only if N(<& — 0) = and + 0) > or, 

equivalently, if and only if S(<E — 0) = — oo and S(<E + 0) ^ 0. 

Proposition 2.9. Alternatively, the ground state energy is given by the zero entropy: 

<B Uii (A,n) =£? w>|t (A,n,0). 

The latter characterization is essentially due to our choice of the inverse function E(S) given 
by (12. 5p and would not be valid for another choice of the inverse, whereas the Proposition 12.81 is 
universal with respect to the particular choice of the function E(S). 
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2.2. Thermodynamic Limit. In this section we discuss the notion of thermodynamic limit, fol- 
lowing the approach of |Rue99j . For sake of completeness and the ease of reading, we repeat here 
the basic definitions related to the notion of thermodynamic limit that can be found in various 
monographs and articles such as |Rue99} ILSSY05| IGri65[ ILL72] . 

First of all, we give a precise meaning to the notion of a sequence of domains tending to infinity. 

Definition 2.10. Let diam(A) be the diameter of A and d^A be the /i-neighborhood of <9A, i.e., 

d h A = dA + B(0,h), 

where B(0, h) is the open ball of center and radius h. 

Definition 2.11. The sets A tend to infinity in the sense of Fisher if 

lim |A| = +oo 

and there exists a "shape function" tc such that 

lim Ti (a) = 

and for sufficiently small a and all A 

|<9 ad iam(A)A|/ |A| ^ ir(a). 

In what follows, we will always assume that A — > oo in the sense of Fisher. 

Remark 2.12. Consider a sequence of rectangular domains. The fact that they tend to infinity in 
the sense of Fisher is equivalent to say that all their sides tend to infinity at a comparable speed, 
i.e., 

d 

J^[[0, Lj] ->■ oo 



i=i 



minj Lj — > oo, 

1 ^ minj Lj/ max, Lj > 1/C. 



Definition 2.13. The limit A — > oo, n/\A\ — > p, where p is a positive constant (density of particles), 
is called the thermodynamic limit. 

Usually one is interested in extensive quantities per particle or per unit of volume (that is the same 
thing up to a multiplicative constant due to Definition 12.131) while considering the thermodynamic 
limit. 

Definition 2.14. Let X W (A, n; V) be a random variable that depends on a domain A, a number of 
particles n and on a set of parameters V. We say that X U (A, n, V) admits the thermodynamic limit 
if the limit 

lim X ^ K ^ V ) 

A->oo n 

n/|A|-T>p 

sm 

exists in some sense with respect to randomness uo (almost sure, in probability, in L 2 ). Here £[P] is 
a certain limiting procedure for the parameters V, i.e., it determines the way how the parameters 
V evolve when A and n go to infinity in the thermodynamic limit. For example, see (13. 3p . where 
an extra parameter is entropy S, and the limiting procedure for the entropy reads as it should tend 
to infinity linearly with the number of particles and/or the volume of the system. 
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In thermodynamics, some commonly used quantities (such as internal energy, for example) are 
assumed to be extensive, i.e., additive with respect to volume. The existence of the thermodynamic 
limit is the mathematically rigorous way of verifying the above assumption. Thus, it is one of the 
fundamental questions of statistical physics. Some authors go even further and refer to the question 
of existence of thermodynamic limit purely as "existence of thermodynamics" |LL72] . 

In what follows, we will be primarily concerned with the existence of the thermodynamic limit 
for the energy E u (A,n, S) with S/n — » u ^ and, in particular, the ground state energy <B u (A,n), 
i.e., for a = 0. 

3. Main Results 
Throughout this section we work with Dirichlet boundary conditions 

H u (A,n) = H5(A,n) 

and we omit the explicit indication D in notations. We give a series of statements concerning 
the existence of the thermodynamic limit for the model of interacting quantum particles in random 
media, which was introduced in Section I2TT1 Basic properties of the thus defined limits are discussed. 
We shall need some assumptions on the model that we introduce now. 

Pair translation invariant interactions. The interactions are by pairs and are invariant 
under translations if for all n G N 

W„(x) = Yl Uix'-x^, (PI) 

where U is a function on V. We also assume that pair interactions are symmetric: U(x) = U(—x), 
xeV. 

Tempered interactions. Assume (PI) and that there exist Ro > 0, A and A > d such that for 
all \x\ ^ Ro 

\U(x)\ ^ A\x\- X . (PTI) 
This condition (together with an additional assumption that U is integrable in a neighborhood of 
zero) guarantees that interactions are of short range, i.e., 

\U(x)\dx < +oo. 

The temperedness or similar conditions on the behavior of the interactions at the infinity have 
been used by various authors such as Leon van Hove, Joel L. Lebowitz, Robert B. Griffits and, in 
particular, Michael E. Fisher and David Ruelle. The reader is referred to |Fis64] . |Rue99] , [FR66J. 
|Leb76j . [Gn65] . 

Remark 3.1. The above assumption of temperedness of interactions can be physically motivated 
by the following argument. Consider electrons in metal or semiconductor as a reference system. 
Though electrons interact via Coulomb potential (~ 1/r) in vacuum, the situation is different in 
metal where each electron is surrounded by a "cloud" of other electrons and lives in a grid of ions. 
This leads to what is called screening of Coulomb potential in metal (see [AM76, Zag98|) and results 
to the effective interaction potential of the form 

U(r) = —exp(-r/X). (3.1) 
r 

The interaction is between quasiparticles "electron+cloud", that are called plasmonsj^ 
1 The potential (|3.ip is called Yukawa potential, though it usually arises in a context of nuclear physics. 
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Lower-bounded one particle Hamiltonian. The one-particle random operator is bounded 
from below uniformly with respect to randomness u: 

3C > 0, such that H w ^ -C for all uj G to. (LB) 

Notation 3.2. We write N n = {1, . . . , n}. For an index set / = . . . , i n } C N we write 

x 1 = (x h , . . . , x in ) G R nd 

for the vector of the coordinates of the particles enumerated by /, where the elements are ordered 
in a nondecreasing fashion: i v < i q if p < q. 

Definition 3.3. Let Ji U Ii = N ni+ „ 2 , \Ij\ = rij, be a partition of n\ + n 2 particles in two disjoint 
subsets. The term of interaction between the particles I\ and I 2 is given by 

W hj2 {x n ^) = W ni+n2 {x n ^) - W ni (x h ) - W n2 (x h ). 

Repulsive interactions. The interactions are repulsive, if for all I\, I2 as in Definition 13.31 it 
holds 

W Iuh > 0. (Rep) 
If one assumes (Rep) and that there are no self-interactions: W\ = 0, then for all n G N 

W n {x N ") ^ Wn^a*"- 1 ) + W x (x n ) > Wn.^x^- 1 ) > . . . > 0, x Nn G R nd . (3.2) 
If one also assumes (PI), then (Rep) is equivalent to say that 

U ^ 0. 

Stable interactions. The interactions are stable if there exists B > 0, such that for all n G N 

W n (x) > -nB. (SI) 

By (13. 2p . repulsive interactions are stable with B = 0. The stability of interactions for various 
models is widely discussed, in particular, in [FR66J. 

Compactly supported interactions. Using the notations of Definition 13. 3[ the interactions 
W have compact support if there exists Ro > such that 

W Iuh (x n "^) = (Comp) 

for all x Nn i+ n i G V ni+n ' 2 such that dist(x 71 , x h ) ^ #0- 

Remark 3.4. Obviously, for pair interactions, compact support is stronger than temperedness, i.e., 

(PI) + (Comp) (PTI) with A = 0. 

Let us now discuss the physical validity of the above assumptions. For more details on classical 
electrodynamics, see, for example |Jac75] and for the electrodynamics of continuous media, see, for 
example [LL60J. 

• The model of pair translation invariant (PI) repulsive (Rep) interactions is natural for a 
description of identical quantum particles such as electrons. 

• The condition of temperedness (PTI) might seem more restrictive at first glance, but is 
usually circumvented as described in Remark 13. II by replacing actual interactions by screened 
interactions and bare electrons by quasiparticles. 

• The condition of compactly supported interactions (Comp) is a technical one and allows 
us to treat interaction of higher order than pair (triple, etc.). However, even short range 
Yukawa interactions (13. ip are not compactly supported. 
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• The repulsive nature of interactions between identical particles (Rep) is widely accepted. 
Though, mathematically only the condition of stability (SI) is needed. Further discussion 
of stability condition and examples of catastrophic, i.e., not stable, potentials may be found 
in |Rue99j . 

• Finally, the lower boundedness of the one-particle operator (LB) seems a natural basic 
assumption. 

The following theorem is the main result of this paper on the existence of thermodynamics for 
the model described in Section [2j 

Theorem 3.5 (existence of thermodynamic limit). Suppose that the one particle operator is lower 
bounded (LB) and that the interactions are stable (SI). Let also any of the following two cases hold: 

(a) interactions are translation invariant and by pairs, i.e., they satisfy (PTI) 

(b) interactions are compactly supported, i.e., they satisfy (Comp). 
Then, the energy per particle admits thermodynamic limit, namely 

K ' ' ' ^£(p,a) as A ^oo, — ->• p, ► a, 3.3 

n \A\ n 

where p > and o ^ 0. The convergence takes place in L 2 (Q) in case (ja| and in L 1 (f2) and co-almost 
sure in case (jb]). The limiting energy density £(p,a) is defined by (13. 3p . is a non-random function 
( does not depend on u) and the limit is the same if both conditions (jsj) and (jb]) are satisfied. 

The energy density has the following basic properties. 

Proposition 3.6 (critical density of particles). There exists a critical density p c G [0, +oo] such 
that 

£(p,a) < +oo, ifp<p c , 
jE{p,<r) = +oo, ifp>p c , 
for all a ^ 0. 

Proposition 3.7 (energy density properties). The energy density S(p,a) is 

(a) a convex function of variables (p -1 ,^); 

(b) a nondecreasing function of p and o; 

(c) a continuous function in the region {0 < p < p c } x {a ^ 0}. 

Corollary 3.8. The energy density S(p, a) admits an inverse cr(p, £). The latter is convex upwards 
with respect to (p _1 ,£ ) and is nondecreasing in 8 for any fixed p. p 

Next we state a reciprocal result exchanging the roles of energy and entropy (the proof follows 
|Gri65j ). 



Theorem 3.9 (existence of thermodynamic limit for entropy). Let the conditions of Theorem^ 
be satisfied. Then for < p < p c and S G Ran£(p, •) 

S u (E,A,n) , . n E 

>a{p,t) as A — >■ oo, — — ■ — > p, \t. 

n \A\ n 

The convergence takes place in the same sense as given by Theorem \S.5\ 

Remark 3.10. The condition that the energy belongs to the image of the function S(p,-) is crucial. 
One might remark as well that due to monotonicity and convexity properties of £ , either £(p,-) 
const identically, or Ran£(p, •) = [inf £(p, •), +oo). 
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4. Proofs 



This section is mainly devoted to the proof of L 2 -convergence (case (jsj) of Theorem I3.5p . The 
basic ideas were inspired by |Rue99] and |Gri65j. though the crucial difference is that instead of 
translation invariance of one particle operator (which is free Laplacian for both of the above works) 
we have ergodicity, i.e., covariance with respect to a family of measure preserving transformations 
of the probability space. 

We assume (LB), (PTI) and (SI) throughout this section, except for Subsection 14.61 where 
different assumptions will be made. We also recall that the Dirichlet boundary conditions are used, 



4.1. Subadditive Inequalities. Subadditive inequalities play the key role in our proofs. The basic 
idea for all the proofs for existence theorems in this paper (and many others: see, for example, 
|Rue99} ILL72j. IGri65] ) may be summarized as: 

• find a subadditive type inequality, 

• use the existing or prove an analog of subadditive ergodic theorem that guarantees the 
convergence. 

Next is the core lemma that gives the subadditivity of energy. 

Lemma 4.1 (Test function construction). Suppose (IAD) and (PTI) are satisfied. Let the statis- 
tics Jj G {0, +, — } be fixed. Consider domains A 1; A 2 such that dist(Ai, A 2 ) ^ r ^ Rq and functions 
(fj G -fjj 3 (Aj), j = 1,2, with energies below Ef. 



Then, using tp\, ip 2 , one can construct explicitly £ G n2 (Ai U A 2 ), a function ofni + n 2 particles 
defined of a unified box A\ U A 2 with energy below E\ + E 2 + An\n 2 r~ x : 



Remark 4.2. The construction of a test function is explicit in the proof that follows. 

Proof of the Lemma. We consider the extensions of the functions ipj, j = 1, 2, by zero on (A!UA 2 ) nj , 
which we also denote by (fj. Remark that (14. ip implicitly contains the fact that (fj are zeros on the 
respective domains boundaries (due to Dirichlet condition) so that the zero extension is a natural 
operation. These extensions obviously preserve (anti) symmetry when (j G {+,—}. Consequently, 
one has ipj G ^^(Ai U A 2 ) for any initial choice of (j. 
We study each statistics separately now. 

Boltzmann statistics: Take 



i.e., H L 




{H^k^nfitpjitpj) ^ EjWipjW 2 , j = 1,2. 



(4.1) 



(^ iS (A 1 UA 2 ,n 1 + n 2 )C,C) ^ (E 1 + E 2 + An 1 n 2 r- X )\\(\\ 2 . 
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Then, using (PTI) 

<^(A 1 UA 2 ,n 1 + n 2 )C,C> 

= ^{^(A 1 ,n 1 )®l n2 + l ni ®iJ w (A 2 ,n 2 ) 

+ W inter {K x , m; A 2 , n 2 )}i fl ® <p 2 , Vl ® <p 2 
< (iy w (Ai,ni)^i,(^i) \\ip 2 \\ 2 + (#u,(A 2 ,n 2 )y? 2 ,v9 2 ) ||^i|| 2 



/ ^ \ <42) 

+ ( U{x % - x J )v?i ® y? 2 , fi <2) y? 2 

\ i=l...ni 
J=n.i+l...ni+n 2 

^ (E 1 + J E; 2 + An 1 n 2 r- A )||C|| 2 , 

where W mter is the potential of interaction between n! particles in Ai and n 2 particles in 
A 2 (see Definition 13. 3p . In the last inequality, we used the temperedness of pair interaction 
potential U. As in each of riin 2 terms, x l G A 1; x J G A 2 and dist(A 1; A 2 ) ^ r ^ R Q , then by 
(PTI) 

dxf/(x i -^')|C(x)| 2 ^Ar- A / / dx|C(x)| 2 ^ Ar- x \\C\\ 2 , 

A™ 1 A™ 2 

where we used the shorthand notation x = x N "i +n 2 for the vector of all particles' coordinates. 
Bosons: We construct 

C(x N «i + " 2 ) = J2 fiix 1 )^^^ 1 ). (4.3) 

card I=n\ 

This function is symmetric with respect to coordinate permutations. Moreover, the terms in 
the sum (14. 3 p are mutually orthogonal in L 2 ((A x U A 2 ) ni+ri2 ); hence, 

ni + n 2 \ 2 2 

Vill \m\\ , 



ni 

and the operator H w , + (Ai U A 2 ,ni + n 2 ) preserves this orthogonality: 

ifw.+CAi U A 2 ,m + n 2 ) • ^{x 1 )^^^ \ J ) 1 ^i(ar J )^ 2 (x NB1+naXJ ), I ^ J. 
Consequently, rewriting (14. 2 p for £ given by (14.31) and using the orthogonality, one obtains: 

(i? W)+ (A 1 + A 2) n 1 + n 2 )C,C) < ■ (#i + ^2 + Anm 2 r- A ) ■ ||^|| 2 1| p 2 || 2 

= (E 1 + E 2 + An 1 n 2 r- X )\\(\\ 2 . 
Fermions: The construction for fermions is similar to that for bosons. We define 

^ x N„ 1+ n a ) = (-l)^ li -<pi(x I )<p 2 (x^+^\ I ). (4.4) 

/CN ni+n2 
card I=n\ 

This function is antisymmetric with respect to coordinate permutations. The remaining part 
of the proof follows exactly that for the bosons. 

□ 

Remark 4.3. The construction (14. 4p is a generalization of the Slater determinant |Gre07] . 
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Remark 4.4. The Dirichlet boundary conditions are crucial for the proof as they provide a zero 
cost (canonical) extension of functions from Domif w (Aj, rij) to a larger domain Ai U A2 without 
changing the norm. At this moment, we are not able to prove an analog of Theorem 13.51 (essentially 
we need an analog of Lemma 14. ip for Neumann or periodic boundary conditions. 

From now on, we omit the statistics sign jj in the notations. Everything that follows is valid for all 
the statistics. However, one should be warned that quantities (such as limiting values) may depend 
on the statistics. 

Proposition 4.5. Let the interactions W be tempered (PTI). // dist(A 1; A 2 ) ^ r ^ Rq, then 

iV u (Ai U A 2 , m + n 2 , E l + E 2 + A ni n 2 r- X ) ^ N u (A u n 1} E 1 )N UJtE (A 2 , n 2 , E 2 ), (4.5) 

S ul (A 1 UA 2 ,n l + n 2} E 1 + E 2 + An l n 2 r- x ) ^ S W (A X , n u E t ) + S W (A 2) n 2 , E 2 ). (4.6) 

Proof. The proof of (14. 5 \ is done using the variational principle for eigenvalues of i^j(A, n) and the 
function £ from Lemma [4. II as a test function. Taking the logarithm, one obtains (14. 6p . □ 

Proposition 4.6. Let the interactions W be tempered (PTI). 

(a) Take S\, S 2 such that exp Sj G N, i = 1,2. If dist(A 1; A 2 ) ^ r ^ Rq, then 

E u (Ai U A 2 , m + n 2 , Si + S 2 ) ^ S w (Ai, ni, Si) + E U {A 2} n 2 , S 2 ) + An x n 2 r- X . (4.7) 

(b) Tafce Si such that exp S, G N, i = 1, . . . , m, and domains Ai, . . . , A m at mutual distances 
greater than r ^ Rq. Then 

tm m m \ m /t / m \ ^ 

|jA i ,^n i ,X;S i U^^(A i ,n i ,S i ) + - I ) r~ A . (4.8) 

t=l 1=1 t=l / i=l \i=l / 

Proof. The inequality (14. 8|) is an immediate consequence of (14. 7J . The latter is obtained by taking 
E u) (A l U A 2 , m + n 2 , •) of (USD and using (gSJ. □ 

4.2. L 2 -convergence on a Special Sequence of Cubes. In this section, we will construct a 
special sequence of cubes Ajy in configuration space V, on which the existence of thermodynamic 
limit will be proven. The idea is inspired by |Rue99]. 
Let 9 be a number that satisfies 

1 < 2 d/x < 9 < 2 

and let 

~ L > R = T^> <«> 

where 5 > is a constant that will be fixed later. For an integer N ^ put 



(4.10) 



so that Ljv £ 2Z, and define the cube A^ by 

A N =[-L N /2,L N /2] d CV. 

Remark that the vertices of A^r are at integer points. According to ( 14. 10p it is possible to place 2 d 
translates of A^v (cubes A$) inside Ajv +1 at mutual distances at least 

R N = L N+1 - 2L N = 6 N {2 -6)R + e = 6 N (R + 5) + e ^ R , 
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where e G [—4, 2] is the error due to the rounding procedure. The constant 5 is chosen to compensate 
a possibly negative error term e, so that the last inequality holds true. It suffices, for example, to 
choose 5 = 4. 

We remark that cubes are explicitly given by 

A («) _ A ,AN) 
A N - Ajv + 7j , 

where 

tP° = Ln+1 ~ Ln ■ « e % d , ei = (±l,...,±l)GM rf , i = l,.. .,2 d . 

Remark 4.7. It is important that the translation vectors 7^ are integer, because it ensures that 
the restrictions of the random potential V u to A^ for different i are connected by the covariance 



relation (ETIjl . 

The function i^(A,n, 5) satisfies the following monotonicity properties. 

Lemma 4.8. For fixed u andn, the energy E u (A,n,S) is 

(a) a nondecreasing function of S, 

(b) a nonincreasing function of A. 

By Lemma 14.81 and the almost-subadditivity condition (14. 8 j) we obtain for 5; in log N that 

(2 d 2 d \ / 2 d 2 d 2 d 

i=i i=i / y=i i=i i=i 

/ j \ 2 



i=l \ i=l 



A 

N 



2 d A / 2d Y 

\(N) M ( A ^ n *> + -77 Z ni I R 

i=l li \ i=l 



A 

N 1 



where {T 7 } 7gZ d is the family of ergodic transformations of Q, that were introduced in Notation 12.21 
In particular, for S G log N 



A 



E u (A N+1 ,2 d n,2 d S) (mM( A N,n,S) + -{2 d n) 2 R- N \ (4.11) 

1=1 ^ 2 

Let now p and o be positive numbers such that 2 N ° d pL d and exp ^2 Nod apL d j are integer, for a 
sufficiently large integer Nq. Plug in (14. lip 

n N = 2 Nd pL d , S N = an N = 2 Nd paL d (4.12) 

for N ^ No. Remark that 

n N /\A N \ p, S N /n N = a, 2 Nd L d /\A N \ -» 1 as iV -» +00. 
We introduce the following sequence of random variables 

^(w) = 2~ Nd (E U (A N , n N , S N ) + (B + C)n N ) , (4.13) 



14 NIKOLAJ A. VENIAMINOV 

where B is the constant from (SI) and C is the constant from (LB). By (14. lip this sequence satisfies 
the inequality 

X N+1 (u) ^ 2- d J2X N (T nv)M) + G N 
i=i 

with 

G N = Ap 2 L 2d 2^ N ^ d ' l R N \ 
In order to show the convergence of the sequence Xn{uj), we establish the following proposition. 

Proposition 4.9. Let Xn be a sequence of nonnegative random variables on a probability space 
(Q,F), Xn(u) ^ 0, such that for each N there exists a family of probability preserving transforma- 
tions of Q, T' ■ , i G In, with card-Zjv < +oo, such that the variables Xn ° t\ , i G In o,re i.i.d. 
(independent identically distributed). If the sequence X N satisfies 

X n +i(uj) ^ — — V X n (t^ N) uj) + G N , where V \G N \ < +oo, (4.14) 

card/jv *rf 

«e/jv n=i 

then there exists a constant X (that does not depend on u) such that 

X N L2 > X. (4.15) 

Proof. Since the terms at the r.h.s. of (14.141) are identically distributed, after taking the expectation 
one obtains: 

EX N +i ^ EX N + G N - (4.16) 

Consider the sequence 

JV-l 

c N = EX N - Gi. 

i=i 

Obviously, (I4.16P guarantees that Cn+i ^ Cjy. Consequently, this sequence converges: Cn > Cq. 

AT— >oo 

Thus, as soon as the sum J^Ja Gi also converges, EXjv admits the limit that we denote by X: 

EX N > X. 

N->+oo 

Consider now the variance of X N , E(X N - EX N ) 2 = E(X N ) - (EX N ) 2 . We will show that this 
variance tends to zero as N goes to infinity. By (I4.14p 

2 



v[T^ n) uj) + G N 



(card/jv) 2 £f V 1 ) (cardi^ 



l£l N 
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Taking the expectation and using the fact that ( X^(t^ ) ) are i.i.d. for fixed N, we find 

i¥=3 



E + <^ E + ^ E + « 

1 ^v2n , (card /*) (card , 2 n _ ^ 2 



E(4) + ^ t^tt^T HEX^ 2 + 2GW ■ KX N + G N . 



card^Af (card In 



So, we get 



< —JtHXn) + ( 1 ?n- J ( EX ^v) 2 + 2G V • EX,v + G^. (4.17) 

card j at \ card i at / 

By the Schwarz inequality 

(EX N ) 2 <: E(X 2 N ). (4.18) 

Hence, using (14. 17ft . we obtain 



E{X 2 N+1 ) <: E(X N ) + 2G n ^/e(X 2 n ) + G 2 N = ^E(X 2 N ) + G r 



Finally, 



e(x 2 n+1 )<:Je(x 2 n ) + g n . 



Arguing as for the expectation, this implies that the sequence E(X N ) converges. Taking the limit 
in (I4.17P and using the fact that Gn — > as N — > +oo, we obtain 

lim E(X 2 N )<:X 2 , 

but the Schwarz inequality (14.181) immediately gives the reciprocal estimate. We conclude that 

lim E(Xl) = X 2 

and the variance E(X^ — EX^) 2 tends to as N goes to infinity. This proves f)4.15p . □ 

As an immediate consequence of this general statement we obtain the existence of the thermo- 
dynamic limit on the special sequence of cubes Ajy. 

Corollary 4.10. Suppose (LB) ; (PTI) and (SI). Then the thermodynamic limit for the energy 
E u on the sequence of cubes in the sense of L 2 ^ exists, i.e., 

E 0J (A N , n N , S N ) l 2 c , \ 

where £(p,cr) is defined by this limit and is called the limiting energy per particle or the energy 
density. 
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Proof. We only need to prove that the random variable Xn{oj) introduced in (I4.13P is nonnegative. 
Recall that 

Hu{k N ,n N ) = ^H$(k N ,n N )+W nN . 

i=l 

Moreover (see also (12. 2p ). 

H®{A N ,n N )^-C, i = l,...,n N , 
by the lower boundedness of the one particle Hamiltonian (LB), and 

W nN ^ -Bn N 

by the stability of interactions (SI). Thus, 

H u (A N ,n N ) > -(B + C)n N , 

and consequently 

E U {A N , n N , S N ) + {B + C)n N > 0. 

□ 

4.3. Critical Density of Particles. We now discuss the finiteness of the thermodynamic limit 
that was announced in Proposition 13.61 

Being essentially attained via a nonincreasing sequence, the limit is finite if and only if there 

are finite terms in the sequence — — — — — — — — . In other words, if for a sufficiently big N, 

n N 

E LU (A N , n N , Sn) is finite, then S(p, a) < +oo. 

The situation when the operator H^{A^,n^) doesn't possess an increasing sequence of eigenvalues 
may arise, according to variational principle, if and only if a subspace of functions ip such that 
{H u {A^,riK[)(p,(p) < +oo is of finite dimension. But this last condition is possible only if the 
interaction potential W n takes the value +oo and if there are too few configurations with a finite 
interaction term, i.e., 

mes&{{x x ,...,x nN ) eA n N N | W n (x\ . . . , x nN ) < +00} =0. (4.19) 

As a model case, suppose that the interactions are by pairs (PI) and that the pair potential U 
represents hard cores of radius r (see |Rue99j ): 

U(x 



= +00, |x| ^ r , 
7^ +00, \x\ > r . 



In this case, the condition (I4.19P is satisfied if there isn't enough space for balls of radius Tq/2 
with centers in the domain Ajy. In other words, define the set of denied spacial configurations of n 
particles by 

S" = {(x 1 , ...,x n ) e M. nd , such that \x % - x j \ < r for some i ^ j} . 
Then the Hamiltonian is defined on £«(V n \ S™ ), instead of Ljj(V") and it may happen that 

meas (V n \ = 0. 

The last observation suggests that there exists a critical density of particles p c such that the energy 
density S(p, a) is finite for p < p c and infinite for p > p c . For example, for the case of hard cores, 
this is the closed packing density. Note that p c = +00 if the interaction potential takes only finite 
values. 
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4.4. Properties of the Energy Density. Before proceeding with the proof of the existence of 
thermodynamic limit for general domains, we establish some properties of the energy density £ (in 
particular, we prove Proposition 13. 7p . 

From now on, we assume that p < p c . Note that till now the function £(p,a) is defined (as a 
limit of a sequence) only for particle densities p and entropy densities a of the form 

rai \ogm 2 
P = 2W a = (420) 

where mi, mi and Nq are positive integers. In all the statements that follow we implicitly assume 
that p and a satisfy (14.20)) . 

Proof of Proposition \3. 71 ([a]) . The convexity of the limiting function is an immediate consequence 
of the almost-subadditivity (14.71) . □ 

Proof of Proposition \3. 71 ()b]) . The monotonicity is given by Lemma 14.81 □ 

Proposition 4.11. The energy density £ is locally bounded on the plane of parameters < p < p c 
and a ^ 0. 

Proof. Let < p\ < p 2 < p c and o~q > be of the form ( 14.20)) . We shall show that £ is bounded 
in the region A = {px ^ p ^ p 2 } x {0 ^ a ^ o"o}. First remark that the number of particles 
hn = 2 Nd L d p with p ^ p 2 can be represented as 



2 



Nd 



n N = J2 n o j) > ( 4 - 21 ) 

where G j L d p , L d p + 1 j. Obviously, such a representation depends on p. On the other 

hand, the bound ^ L d p 2 +1 =: n™ ax depends only on p 2 . This representation can be obtained 
as the result of a consecutive division of the domain A^v in sub-domains (each time we divide the 
domain in 2 d parts) until one obtains the domains Aq \ In (14. 9)) . choose L sufficiently large so that 
L d p ^ 1, i.e., there is at least one particle in each sub-domain Aq^. 
Let us denote by S* the smallest number belonging to logN that is larger than S: 

S* = inf{Q ^ S, exp Q e N}. (4.22) 

Then one calculates: 

A 
2 

j=l m=0 

Since S N /2 Nd ^ L d p 2 a and sC n^ ax , we can deduce that E u (f$\r$\ (S N /2 Nd )*^ is bounded 

uniformly with respect to N and u>. The proof is done by a trivial bounding of the potential (as the 
number of terms in the potential is bounded) and by the application of Weyl asymptotic. Dividing 
P3| by 2 Nd , we finish the proof. □ 

Proof of Proposition ^. 71 (jcj). Having established Proposition 14.11) it is sufficient to apply a standard 
argument due to Jensen (see, for example, |PS98j). □ 



E U {A N , n N , S N ) ^J2 E » ( A o J) > n o\ (S N /2 Nd Y) + -n N nZ x £ 2^ d R~ 2 . (4.23) 
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Corollary 4.12. Suppose p < p c is of the form ( I4.20p . o\ and 02 are fixed. Then the L? -convergence 

n- N l E„ (A N , 2 Nd L d p, 2 Nd L d pa) -> £{p, a) 

is uniform in o G [01,02] because the pointwise convergence of monotone functions to a continuous 
function on a compact interval implies the uniform convergence (Dini's theorem). 

From this corollary we deduce the following proposition, which weakens the restrictions on the 
way the entropy must go to infinity in the thermodynamic limit. Instead of a very specifically 
chosen sequence Sn as in Corollary I4.10[ we need only the linear dependence of entropy on the 
number of particles. 

Proposition 4.13. Let p > and a ^ 0. Analogously to (14.1 2ft . construct a sequence of integers 
Un = 2 Nd L d p. Let also Sn be a sequence such that > a. Then 

n N E u (A N ,n N ,S N ) > £(p,o~). 

N— >+oo 

4.5. Proof of Theorem 13.51 ^aj) (L 2 -Convergence for General Domains). Now we are ready 
to show the existence of the thermodynamic L 2 -limit in full generality. First of all, we will establish 
that 

limsup Eaj(A,n ' ^ ^ S(p,a). (4.24) 

A— >oo n 

Let po > P of the form (I4.20P be close to p: 

^ = ^k' K,iVo)GN 2 . (4.25) 
2 N ° d L d 

The representation (I4.25P is not unique. Among all the representations for a fixed po there exists 
the minimal one, i.e., where Nq is minimal: po = n m i n 2~ Nmmd L~ d . We write: 

n = mn + r , ^ r < n . 
For a fixed po we choose Nq in (I4.25P as a function of n such that 

L d N Jn^0, L x N Jn-^ +00. (4.26) 
By Definition 12.111 a sufficiently large A contains m disjoint cubes with sides £Ln , where 

l/d 

i 1 >l\ \ 



Consequently, A contains translated cubes A^, i = 1, . . . , m, at mutual distances at least (£ — 1)Ln - 
By one has 

m-l . 

E U (A, n,S)^Yl E ^ A n > n ^ S K m - X )) + E ^N^n + r , 0) + -n 2 (£ - 1)" A L^. (4.27) 



i=l 

We treat now the term E LJ (A^,n + r ,0). By the domain division procedure similar to that 
described in Section 14.21 we can reduce A^ to the union 

of 2 ( N "-N min )d = no / nmin (this is an 

integer) translates of A^ min . We obtain 

^(A^no + rcOX £ M^ min , n mm + r j} 0) + -(2n ) £ 2^ d R~ m 2 (4.28) 
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where ^ r,- < n min and j G {1, . . . , n /n min }. Clearly, 

^>(AivL„> n min + r j ; 0) < C\n 2 min 
where the constant C\ is uniform in j and u. By (I4.28P one obtains 

Finally, from (I4.27P one deduces that 

rT^A, n, 5) < ^ -^(41 "o, S-/(m - 1)) + C 3 ^^ + - 1)-*-=^. 

Note that uq/u ^ 1/m — » and n/L^ o — >■ according to (I4.26p . Thus, 

lim sup n" E W (A, n, S) ^ £ (po, cr) • 

A— >oo 

Approaching p from above by p , one gets ( I4.24p . 
The proof that 

lim inf n~ l E u (A, n, S) ^ £(p, a) 

A— i>00 

is done exactly as in |Rue99] . pp. 47-48. □ 

4.6. L 1 and Almost Sure Limits. In this section we assume (LB) and (Comp). We show 
that if the interactions are compactly supported, the convergence to the thermodynamic limit can 
be improved. The proof follows that of the L 2 -convergence, so we only indicate the necessary 
modifications. 

Let us introduce the following random variable. For all n G N, S G R and domain A C M. d we set 

f u (A,n,S) = E u (A,n,S*), 

where A is the .Ro/2-interior of A, i.e., 

A = {x G A, dist(x, dA) > R /2}, 

and recall that S* is defined in (|4.22|) . We make the following observations before giving a subad- 
ditivity condition. 

Lemma 4.14. Let A and B be two domains in V. 

(a) If Ad B, then Ac. B. 

(b) If Af) B = 0, then dist(A, B) ^ R . 
Let also x, y G M. 

(a) If x ^y, then x* ^ y*. 

(b) (x + y)* ^ x* + y*. 

We modify now the subadditive inequality (14.71) . 

Proposition 4.15. Let the interactions W be compactly supported (Comp). Let ni,ri2 G N. If A 
and B are two disjoint domains in V, An B = 0, then 

UA UB,m + n 2 , S ± + S 2 ) ^ U(A, n u St) + f u (B, n 2 , S 2 ) . (4.29) 

Proof. As the interactions are compactly supported, we get (14.71) by the same manner as in Propo- 
sition H]6] but without the interaction term (with A = 0). It remains to use Lemma [4.141 and the 
monotonicity of energy with respect to entropy in order to get (14. 29 p . □ 
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Thanks to the subadditivity (14.291) . we prove the convergence in L 1 and the almost sure conver- 
gence. 

Proof of Theorem \3.5\ (jb]). In order to prove this type of convergence, which is stronger than that in 
the part (jsj) of the theorem, it is sufficient to modify Section 14.21 Everything what follows remains 
true without any modifications. 

In Section I4T21 we change the definition of cubes by taking 



L 



N 



2 N L 



Rn-S 



in a place of ()4.10|) . where S is a fixed positive constant. This guarantees that one may put exactly 
2 d translates of in a cube A^r+i at distances at least Ro for a properly chosen 5. The lower 
boundedness is given by (LB) and (SI): 

UA,n, S) 



>-B-C. 



71 



Next, we apply the multidimensional subadditive ergodic theorem (see, for example, |Smy76| ) and 
obtain the L 1 - and almost sure convergence of the sequence. □ 



5. Free Particles 

As a complement, we study the thermodynamic limit for the energy density E(p,a) in the case 
of free (noninteracting) particles: 

W = 0. (5.1) 

We remark that the background potential V w remains present. Interestingly, even in this case the 
results are not as trivial as one could have expected. The obtained thermodynamic limits depend 
on quantum statistics. 

5.1. Maxwell-Boltzmann Particles. For particles without statistics we establish the following 
theorem. 

Theorem 5.1. Suppose that the interactions are absent (15. ip and that the particles are of Maxwell 
- Boltzmann statistics. Let E be the almost sure spectrum of the one-particle Hamiltonian H w (l). 

If 

E = suppdN ; (5.2) 

then 

S(p, a) = inf E 

for all p > and o ^ 0. 

Remark 5.2. The condition (15. 2p is satisfied under rather general assumptions on the random po- 
tential Vu (see, for example, [Ves08j). 

In order to prove Theorem 15.11 we will make use of two following lemmas. We assume that the 
conditions of this theorem are verified in the sequel. 

Lemma 5.3. Let u be such that Spec(i7 U) (l)) = E. Let N e N be fixed. Then 

E N (H U (A,1)) -> inf E, A->cx). 
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Proof. Consider, as usual, Dirichlet boundary conditions. Then for almost any £6 1, 

N A (E)/*N(E) } A^oc, 

where N A is the pre-limit density of states, i.e., the counting function of the operator iJ w (A, 1) 
divided by |A|, and N is the density of states of the one-particle operator 
Just by the definition of the counting function 

N 

Na(E n (H lu (A, 1))) = J^yJ — ^ 0' A^oo. 

Moreover, by the monotonicity of the Dirichlet eigenvalues, E^^^A, 1)) decreases as A — > oo and 
thus necessarily converges: 

E N (H W (A,1))\E, A^oo. 
By combining the previous arguments, we find that 

N(E) \ N A (E) < N A (E N (H U (A, 1))) ^ 0, A -> oo, 

which implies 

N(E) = 0. 

Consequently, 

E < inf E. 

Finally, E < inf £ is impossible, because necessarily En(H w (A, 1)) ^ inf S. □ 

The next lemma is a modification of Proposition 14.61 and expresses a subadditive-type condition 
that is even stronger than (14. 8p . It is exactly this lemma that is not valid for bosons and fermions. 

Lemma 5.4. If W = and the particles under consideration are not restricted to any statistics, 
then 

E„(A, m + n 2 , S, + S 2 ) ^ E U {A, n u S,) + E U {A, n 2 , S 2 ). (5.3) 



Proof. The proof follows that of Proposition 14.61 The only modification is the construction of test 
functions in Lemma H~T1 When the interactions are absent, one can place two groups of particles in 
the same box and, consequently, one is not obliged to enlarge the size of a box together with the 
number of particles. □ 

The last idea is not applicable to bosons or fermions, as one cannot guarantee the independence of 
constructed test functions and the orthogonality of terms in (14. 3 p or ( 14. 4 j) is not assured. Moreover, 
a constructed test function may happen to be identically zero for fermions. 

Proof of Theorem \5.1i Because of the subadditivity (15. 3p . the following limit exists: 

3 lim E ^ S) =:C(A,a) 

rwoo n 
A fixed 

in the sense of L l and almost surely with respect to u. 

This is proved exactly in the same manner as the existence of S(p, a) (and in some aspect is even 
simpler). As Dirichlet eigenvalues are monotonous with respect to the domain, the function ( is 
nonin creasing in A and nondecr easing in a (by obvious reasons). Due to the monotonicity in A, we 
find also that for p > : 

£(p,a) ^ lim C(A,a). (5.4) 

A— >oo 
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We remark as well that by (15.31) the function E U} (A, n, an)/n is nonincreasing in n. So one can 
interchange the two limits in the r.h.s. of (15. 4ft to get 

£(p,a) < lim lim — — — - — - — - = lim lim — — — - — - — - < lim lim E W (A, 1, S/n). 

S/n-^cr S/n-^a S/n-^a 

Here (15. 3p was used once more in the last inequality. Finally, by Lemma [5.31 we establish 

lim EjA,l,S/n) -> infS 

A— >oo 

with S and n being fixed. □ 

Theorem 15.11 expresses the fact that the thermodynamic limit for Maxwell - Boltzmann particles 
is trivial in the absence of interactions. Thus, it is indeed the interactions that may possibly render 
the limit being nontrivial. 

5.2. Bosons. We remark that the energy levels for a system of noninteracting bosons coincide with 
energies for Maxwell - Boltzmann particles (see, for example |LL77[ IGre 07j). On the contrary, for 
bosons the combinatorial degeneracy is lifted up by means of the symmetrization procedure (the 
degeneracy due to coincidences like Ei + E$ = E\ + E4 remains) . 

Nevertheless, the ground state energy for bosons is the same as in the previous section, and, 
consequently, the ground state energy per particle converges in the thermodynamic limit to the 
lower edge of the almost sure spectrum of the one-particle operator. 

5.3. Fermions. The situation changes significantly if particles are fermions. For the basic proper- 
ties of a system of noninteracting fermions, we refer the reader once more to |LL77t Ore07j. 

For fermions we know only how to obtain results on the ground state energy S(p, 0). The argu- 
ments we use do not rely on subadditivity properties and, consequently, are valid for any boundary 
conditions. 

The main difference between fermions and bosons is that the ground energy for n noninteracting 
fermions is given by the sum of the first n energies of a one-particle system 

n 

E 1 (A,n) = Y,E k (A,l) (5.5) 

k=l 

and not by n times the one-particle ground energy. The ground state itself is given by the Slater 
determinant 

f2i(A,n) = det C0i(^'))i,j, 

where fa is the eigenfunction of H U (A, 1) corresponding to the energy i£»(A, n, u). 

A comparison with the Laplacian and the use of Weyl asymptotic provide a simple proof that 
the limit S(p, 0) is strictly different from zero for p > if the background potential is nonnegative. 

Proposition 5.5. Suppose ^ 0. Then there exists (3 = /3(d) such that 

£(p,0)>/3p 2 / d . (5.6) 

Proof. By the variational principle, as the potential is positive, we obtain 

^•(A,1)^^(A,1), 
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where E® are eigenvalues of —A. Next, with the help of Weyl asymptotic for the Laplacian, by 
( 15. 5 p we get 

= l -± E k (A, l)> l -± E°(A, (JL) . n Wi > C ^l\ 

fe=i fc=i fc=i Vl 17 11 

which proves (15. 6p . □ 
Remark 5.6. The generalization to the case of lower-bounded random potential is obvious. 

We now compute an explicit expression for the limit 

Um g ,(A,l) + ... + £ „(A,n) 

A— >oo n 

n./|A|— >p 

in terms of the integrated density of states of the one-particle problem. Once more for simplicity, we 
suppose that background potential is nonnegative: V u ^ 0. We need only rather general assumptions 
on the density of states, which we denote by N(E). 

Condition 5.7. The integrated density of states N(E) is a continuous function and defines a 
positive measure dN(E) such that the almost sure spectrum is equal to the support of diV: 

supp diV = E. 

Remark 5.8. The last condition is certainly verified, for example, if the Wegner estimate (W) holds 
for^(A,l). 

Proposition 5.9 (Wegner estimate). Let E be almost sure spectrum of H w . There exists constant 
C > such that, for any Borel subset Icl, 

E(Tr(l/(^(A)))) < C|A| • |/HE|. (W) 

Proof. Wegner estimate is well known for both discrete and continuous Anderson model under the 
assumption that the random variables are i.i.d. and that their distribution is regular |CHK07[ 
IVes08j . □ 

Definition 5.10. Fix a density of particles p. The Fermi energy E p is a solution of the equation 

N(E p ) = p. (5.7) 

Remark 5.11. It may happen that sup E N(E) < p. For example, if one considers a discrete Anderson 
model, then N(E) ^ 1 and for p > 1 the equation (15.71) does not have any solutions. This is due 
to the fact that the density of particles is too big (in other words, there isn't enough space for so 
many particles) to accommodate for n fermions. This situation never arises in a continuous setting. 

Remark 5.12. A solution of the equation ( 15. 7p is not necessarily unique if the integrated density of 
states is flat on the level p. As N(E) is a continuous nondecreasing function, the set of solutions 
is the closed interval [£"™ m , E™ ax ]. From the spectral point of view the open interval (-E™™, E™ ax ) 
doesn't play any role because its intersection with the almost sure spectrum E is empty. 

In this situation we will also use the notation introduced by the Definition I5.10[ meaning E p = 
^jjjmin ^ jjjmax^ we w jj] gee ^ fais convention is consistent with the results. 

Next theorem is the main result of this section. 
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Theorem 5.13. Let p > 0. Then 



S(p,Q) = - / EdN(E). (5.8) 
P Jo 

To give a proof to this theorem, we will need the following crucial lemma, which explains why 
E p given by (15. 7p corresponds exactly to the common physical notion of the Fermi energy^ 



Lemma 5.14. 



E n (K 1, w ) — — >■ E p . (5.9) 

A— »oo 
ro/|A|— ip 



Proof. We denote by the density of states before taking the limit for a one-particle operator: 

N^(E) 



tA/ ^ M(E,H u (A,1)) 



A 



Then by definition, in the thermodynamic limit 



N^E n (A,l,u)) = ^^p. 
On the other hand, by the existence of the integrated density of states we get: 

A— >oo 

We finish the proof by applying the monotonicity argument. Suppose that 



liminf E n (A, 1, to) < E 



mm 
P 



By passing to a subsequence, we find that there exists 5 > such that E n (A, 1, u>) < E r p nm — 5. We 
arrive to a contradiction: 

N ( E rnin^ ^ N^(E n (A, l,w)) < N ul {A{E p ni " - 5)) ->• n - 5) < Ar(E™ n ). 

The last inequality is strict because E p nm is the minimal value of energy such that N(E) = p and 
so for any E above this level the density of states N(E) is strictly smaller. 
Similarly, we show that limsup E n (A, l,u) ^ , so 

£™ n < liminf E n (A, l,u) < limsup E n (A, l,w) < £™ aa; , 

which is equivalent to (15.91) . □ 

Proof of Theorem \5.13l To show ( 15. 8 p we write 

= " = ~ Tr IM A ' !) • l[o,i?„(A,i,-)](^(A, 1))J 

|A| rE n (A,l,uj) i /-Bp 

JU n/|A|->P 

where the convergence is valid because the measure dN A converges weakly to dN, the integration 
limit converges to E p by Lemma 15.141 and the dominated convergence theorem can be applied. □ 



2 Fermi energy is the energy of highest occupied quantum state in a system of fermions at absolute zero temperature. 
Alternatively, for non-interacting fermions, it is the increase in the ground state energy when one particle is added 
to the system. For more details on the concept on the Fermi energy, see [A~M76j . 
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Remark 5.15. The formula (15. 8p admits an alternative form: 



£(p,Q) 



J Ep EdN(E) 
J Ep dN(E) 



which reads as the ground state energy density is the energy averaged from zero to the Fermi energy 
with respect to the density of states. 
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